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Abstract 


This communication presents a model of electrochemical cells developed in order to simulate their electrical, chemical and thermal behavior 
showing the differences when thermal effects are or not considered in the charge—discharge process. The work presented here has been applied to 
the particular case of the Pb,PbSO,4|H2SO, (aq)|PbO2,Pb cell, which forms the basis of the lead-acid batteries so widely used in the automotive 


industry and as traction batteries in electric or hybrid vehicles. 


Each half-cell is considered independently in the model. For each half-cell, in addition to the main electrode reaction, a secondary reaction is 
considered: the hydrogen evolution reaction in the negative electrode and the oxygen evolution reaction in the positive. 

The equilibrium potential is calculated with the Nernst equation, in which the activity coefficients are fitted to an exponential function using 
experimental data. On the other hand, the two main mechanisms that produce the overpotential are considered, that is the activation or charge 


transfer and the diffusion mechanisms. 


First, an isothermal model has been studied in order to show the behavior of the main phenomena. A more complex model has also been studied 
including thermal behavior. This model is very useful in the case of traction batteries in electric and hybrid vehicles where high current intensities 


appear. 


Some simulation results are also presented in order to show the accuracy of the proposed models. 


© 2006 Elsevier B.V. All rights reserved. 
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1. Introduction 


The modeling of electrochemical cells, secondary batteries 
or electrochemical energy storage devices in general has been 
widely studied, there being a considerable number of previously 
published works available [1-20]. The large number of works 
carried out around this topic is partly motivated by the interest 
aroused by electric and hybrid vehicles. For this reason, many of 
these works have been developed as a part of simulation models 
for these vehicles [1—6]. 

These works are sometimes based on empirical relationships, 
and other times on a detailed description of the physical and 
chemical processes that take place in the cell, and even on the 
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development of equivalent circuits [6—8]. In order to develop 
these models various techniques have been used, such as lookup 
tables, lumped parameters or bond graphs. 

Taking the first works published as a starting point, differ- 
ent extensions and modifications of increasing complexity and 
precision have gradually been proposed. In some cases [4,9-11] 
modifications have been made to previously published works 
[3,12], while in other cases, it has been the authors themselves 
[13] who have made suggestions for improving their own models 
[10]. 

There is also a wide variety depending on the goals. There are 
isothermal models, for instance, to describe electrochemical cell 
behavior under normal charge or discharge conditions, models 
that seek to get as close as possible to the limit current condi- 
tions [13], thermal models that include the effect of temperature 
[1,7,8,14,15], and even age on cell performance [1], models 
that study the recombination of gases during the charge pro- 
cess, particularly in lead-acid batteries (VRLA) [10,13,16,17], 
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or models specially developed to describe cell behavior during 
a pulse charge/discharge operation [17,18]. 

A highly interesting alternative for developing these types 
of models is to use the bond graph technique. This technique 
allows for the development of models that may belong to dif- 
ferent physical domains. This is especially useful to integrate 
the models developed into other more complex and multidisci- 
plinary ones as, for example, a complete hybrid vehicle model, 
with all the components like engine, clutch, gear box, chassis, 
and power system, including the battery and the control system. 

In this paper, based on two previous ones by the authors 
[21,26], a model for an electrochemical cell has been devel- 
oped. The model has been applied to the particular case of the 
Pb,PbSO.4|H2SOq (aq)|PbO2,Pb cell. This kind of cell forms the 
basis of lead-acid batteries that are widely used in the automo- 
tive industry as well as traction batteries in electric or hybrid 
vehicles. 

The main goal of this work has been to obtain a model capa- 
ble of reproducing cell behavior for both charge and discharge 
conditions. Therefore, an initial isothermal model has been con- 
structed. In many cases, this isothermal model may be sufficient 
to provide an adequate representation of cell behavior. How- 
ever, in those cases where high charge currents are involved, the 
thermal effects may be considerable. A thermal model has also 
been developed for this case, aimed at providing an adequate 
representation of these effects. 


2. Electrochemistry of the Pb,PbSO4|H,SO4 
(aq)|PbO2,Pb cell and bond graph models 


2.1. Equilibrium potential and open circuit voltage 


Let us consider a cell such as shown in Fig. 1. It consists 
of a lead electrode (Pb) and another of lead dioxide (PbO2) 
submerged in an H2SO0; solution. When the two electrodes come 
together (Fig. la), an electric current flows from one to the other 
due to a potential difference between the two electrodes. The 
reaction giving rise to this electromotive force is: 


PbO2 + Pb + 2H2S04 —> 2PbSO4 + 2H20 (1) 


cathode 


(a) 


galvanic cell 


cath ode_ 


(b) 
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This reaction can be split into two parts. One that takes place 
in the left electrode: 


Pb + S04? > PbSO4+2e7 (2) 
and the other that takes place in the right one: 
PbO? + 4H + S04% + 2e7 —> PbSO4 + 2H2O (3) 


The electrode where oxidation occurs is called the anode 
and the one where reduction occurs is called the cathode. The 
electromotive force of a galvanic cell is defined as: 


E = Es — Ea (4) 
where E, and E, are the reduction potentials of the cathode and 
the anode, respectively. 

If a voltage greater than its electromotive force is applied to 
the cell, the electrochemical process is inverted in the electrodes 
(Fig. 1b) and the direction of the reactions (2) and (3) is also 
inverted. 

The dependency of the electromotive force on the concen- 
tration of substances taking part in the reaction is given by the 
Nernst equation: 


E (5) 


o R T vi 

= E? — — ln (miyi) " 

nF i 
where E° is the normal or standard electromotive force, R 
the molar gas constant, T the electrolyte temperature, F Fara- 
day’s constant, n the number of electrons exchanged, v the 
stoichiometric coefficient, m the molality, and y is the activity 
coefficient. In our case, because of the high concentrations of 
electrolyte used, the solvent activity cannot be considered con- 
stant and equal to the unit and must be included. The values of 
water activity can be obtained, for any concentration of H2SO4, 
from experimental data [22] by adjusting to a polynomial func- 
tion. 

In the case of the electrolytes, owing to the electrically neu- 
tral condition, the concentration of one of the ions cannot be 
varied without changing that of the other, since it is impossible 
to separate the effects of the positive ions from the negative ones. 
This means that the activity coefficients of the individual ions 


cannot be measured experimentally. For this reason the mean 
anode 


electrolytic cell 


Fig. 1. Pb,PbSO4|H2SO4 (aq)|PbO2,Pb cell. 
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ionic activity coefficient y+ is defined as: 


wey xy; v=vt + (6) 


this being the only coefficient that can be measured experimen- 
tally and able to be found in tables. In our case, as the reactions 
of each electrode are to be looked at individually, it is neces- 
sary to have some procedure to estimate the value of the activity 
coefficients of the positive and negative ions separately. Tak- 
ing Debye—Hiickel’s theory for strong electrolytes as a starting 
point, a relation between these coefficients and the mean ionic 
activity coefficient can be established: 


y (v/v+) In y4 
~ 4 /z_-? + _/v4) 
(7) 


where z+, y and z_ are the charges of the positive and negative 
ions (in our case SO4?~ and H* ions) and: 


_ (v/v_)Inye 
 @_ farts) 


In y4 


y = 0); y- =y (8) 


The values of y+ can be obtained for any concentration of 
H2SOa, from experimental data [22] by adjusting to a polyno- 
mial function. 

The value given by the Nernst equation corresponds to the 
reversible electromotive force, that is, when no current is passing 
through the electrodes or when the current can be considered as 
negligible and is called open circuit voltage Under the same con- 
ditions, the potential of each individual reaction in expression 
(4) is known as equilibrium potential. 


2.2. Overpotential 


Under irreversible conditions, that is, when current is flowing 
through the electrodes, the value for potential is different from 
that of equilibrium, the difference being known as overpotential, 
and represented by 7. In our case, of all the different mechanisms 
that can give rise to overpotential, the most important ones are 
those of charge transfer and diffusion, and only these will be 
taken into account. 


2.2.1. Charge transfer overpotential 

Charge transfer overpotential is caused by the charge transfer 
process between the electrode and the electrolyte through the 
double layer. The electrochemical kinetics shows its relation to 
current density by means of Butler—Volmer’s equation: 


i zak z1—a)F (9) 
7 exp RT M exp RT Ni 


where ip is the exchange current density, z the number of elec- 
trons exchanged, a the charge transfer coefficient and 7 is the 
charge transfer overpotential. This relation is shown in Fig. 2 by 
the solid line, while the dotted lines correspond to each of the 
addends of expression (9). 

However, due to the features of the model, current density 
has to be used as an independent variable, and, therefore, the 
inverse relation nų =f! (i) is necessary which, given the form of 
(9), cannot be directly calculated. Nevertheless, as can be seen 


114 


Fig. 2. Current density according to the charge transfer overpotential. 


in Fig. 2, for high values of |7;|, one of the two addends in (9) 
can be neglected, thus obtaining Tafel’s equations: 


em > Fs m>0 
m = —- In- (10) 


e in| > RT, m<O 
RT io 


= oe 


Nt 


For values of /i/ near to iọ and lower, the representation of n 
by Tafel’s equations is inadequate and leads to a considerable 
error in 7. We have chosen to seek a function i=f(7,), with 
behavior similar to that of (9) in the small value zone of /i/ 
and that, at the same time, allows to obtain the inverse relation 


m =f (i). 


2.2.2. Diffusion overpotential 

When current is flowing through an electrode, the so-called 
diffusion layer is created as a result of the transport of the differ- 
ent species taking part in the reaction from or to the electrode. 
This layer is a zone where the concentration gradient is other 
than zero. As a result of this change in concentrations, a differ- 
ence of potential is caused regarding the equilibrium potential. 
The corresponding overpotential produced is known as diffusion 
overpotential. 

The diffusion layer takes on a complex profile and is usually 
substituted by a simpler model, called the Nernst layer, where 
the concentration gradient is taken as constant. By using this 
approximation, the electrochemical kinetics gives the expression 
(12) for diffusion overpotential, according to current density, and 
the density of current limit of each of the species i); diffused 
through the electrode/electrolyte interface: 


RT iN» 
na = — In 1- (12) 
j 


nF ij 
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Fig. 3. Mixed potential formation Ey in the anode. 


where n in this case, is given by: 


n= Sov <i (13) 
j 

and where v; and z; are the stoichiometric coefficient and the 

charge of species j, respectively. 


2.3. Secondary reactions and mixed potential 


Frequently, two or more processes take place simultaneously 
in an electrode, as occurs with the anode where, apart from the 
main reaction, the hydrogen evolution reaction takes place: 


2H* +2e° > Ho (14) 


The current density versus potential behavior curve for the 
electrode, where the two reactions shown take place, is repre- 
sented in Fig. 3. 

The curve on the left shows the charge transfer process in 
the electrode main reaction, where n; is the overpotential of this 
reaction, and Ep (1) is its equilibrium potential. Likewise, 72 is 
the overpotential of the reaction (14) and Eo(2) its equilibrium 
potential. Because the two reactions are independent, the total 
current density i is obtained by adding the densities of the two 
reactions, that is: 


i=i +i (15) 


in such a way, that when i=0, the potential is neither Eo(1) nor 
Eo(2), but an intermediate potential Eyy, which will be called 
mixed potential, resulting in: 


Eo(1) + nı = Eo(2) + m = Em (16) 


Likewise, the following secondary reaction occurs in the 
cathode: 


H20 > 50) + 2Ht + 2e7 (17) 


with oxygen evolution, and with everything commented regard- 
ing the anode remaining equally valid. 


(a) (b) 


Fig. 4. Equilibrium potential representation. 


3. Simulation model 


The simulation model has been developed using block 
diagrams in Simulink, and implementing the equations and 
dynamic model of the battery in this program. In order to under- 
stand the composition of the model, a model with its equivalents 
in electrical components has been developed that represent the 
different mechanisms for accumulating potential energy (com- 
pliances or capacitors), kinetic (inertances or inductions) or 
energy dissipation (resistances or resistors). The structure of 
the model and how the aforementioned electrical equivalents 
represent the phenomena are described in the above section. 


3.1. Isothermal model 


The reversible behavior of the cell in respect of the equilib- 
rium potential has been modeled by using a capacitor (Fig. 4a). 
This element stores free energy, integrating the current inten- 
sity (i), in order to obtain the electric charge (q), and thereby 
obtaining the equilibrium potential by means of Eq. (5), since 
the electrical charge and molality are directly related by Fara- 
day’s laws. Fig. 4b represents the block diagram relative to this 
element. From Fig. 4 onwards, the number or numbers in brack- 
ets placed inside the block refer to the expressions defined in the 
text of this paper. Thus, starting out from a state of initial charge, 
and for a positive current, the charge will gradually increase over 
time, as will the stored electrochemical energy. When the cur- 
rent is negative, the charge will gradually diminish, in the same 
way as the stored energy. When the state of initial charge is once 
again reached, the equilibrium potential will also be the same as 
at the point of departure. 

The irreversible behavior of the cell has been represented 
using a resistance. This element directly relates the current den- 
sity (i) with the overpotential (7) through expressions (9-12), 
as is shown in Fig. 5. As a consequence of the principal of the 
positive production of entropy, the relations between the inten- 
sity and overpotential variables must remain in the Ist and 3rd 
quadrant, as can be seen in Fig. 2. In this way, the potential asso- 
ciated with the resistance will always be positive and represent 
the dissipated electrical potential. 

There should be a capacitor and a resistance for each reac- 
tion. Both elements should be joined so that current intensity is 


ANN- 
i 


a oa 


(a) (b) 


Fig. 5. Overpotential representation. 
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Cy Rs 


Fig. 6. Main and secondary reactions representation. 


the same for both, and the overpotential added to the equilib- 
rium potential. In each electrode, the elements corresponding to 
each of the reactions (main and secondary) should be joined to 
one another so that they conform to expressions (15) and (16). 
Fig. 6 represents the equivalent circuit diagram for the main and 
secondary reaction. 

Finally, the group of elements corresponding to the cathode 
should be joined to those of the anode so that the current traveling 
through both will be the same. The electromotive force of the 
battery will be given by expression (4), where E. and E, are 
the mixed potentials of the cathode and anode, respectively, in 
line with expression (16). Fig. 7 shows the equivalent electrical 
diagram and its representation in a block diagram. 


3.2. Thermal model 


In any chemical reaction, like in those described above, 
besides the electrochemical energy, there are other energy trans- 


Cme Rme 


formations such as thermal energy and mechanical energy. The 
Gibb’s free energy variation will be expressed by: 


dG = -Sev dT + VdP + EdQ (18) 


where dT, dp and dQ are temperature, pressure and charge vari- 
ations, respectively, Srey, V and E the entropy (the part owing 
to the reversible process), volume and the previously defined 
equilibrium potential. 

A special case would arise when the temperature and pressure 
remain constant. It can easily be seen from Eq. (18) that under 
these conditions the changes in the available electrochemical 
energy arise out of the variations in the Gibb’s free energy, the 
reversible behavior of the cell being able to be represented as 
set out in Section 3.1. 

In general, in a cell like the one described here, it can be taken 
that any process goes on at constant pressure. In many cases, par- 
ticularly where there are no very high charge currents, it can be 
taken that the temperature remains equally constant. However, 
when the intensity of the current is high, as happens with traction 
batteries for electric vehicles, the thermal effects can become 
considerable. Under these circumstances, the temperature can- 
not be taken as constant. In this case, the reversible behavior of 
the cell should be represented (maintaining the approximation 
of constant pressure) using a capacitor where, in addition to a 
chemical port, like that used in Section 3.1, should be added a 
thermal port, aimed at taking account of the thermal energy flows 
produced as a result of the reversible part of the heat reaction, 
known as the Peltier effect. 


Cma Rma 


| 


Fig. 7. Full battery representation for the isothermal model. 
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E(0) 


(5)+(19) 
(9+12)+(20) 


Fig. 8. Block diagram of a chemical reaction considering thermal effects. 


S ir 


Via the thermal port is received the temperature value infor- 
mation, on which, in turn, the equilibrium potential depends, and 
which serves to calculate the Srey variable using the following 
constitutive relation: 


dG OE 
Srey = sa) = os). (19) 


A port must be included in the block diagram, in order to 
represent the production of entropy as a consequence of the 
irreversible process (Joule effect), in such a way that the ther- 
mal potential produced coincides with the electrical potential 
dissipated: 


T Sie = nI (20) 


where /is the total electrical current due to the electrode reaction. 
Fig. 8 shows the corresponding block diagram. 

The sum of Srey and Sirr of each of the electrode reactions will 
give rise to heating or, circumstantially, due to the reversible 
process, cooling of the electrolyte. This will act as a thermal 
accumulator and supply the temperature starting out from the 
value of the entropy in line with the relation: 


S 
T = Tef exp (= ) (21) 
V 


where Tref stands for the temperature in the initial conditions of 
the system, m the electrolyte mass and c, the specific heat per 
mass at constant volume. For small variations in temperature, 


i 
| Be (5)+(19) 


| _}(9+12)+(20)] 


and taking the calorific capacity of the electrolyte as constant, it 
can be established [23] that: 


T = Tref = Tref 


(22) 
m cy 
If the temperature Tis different from the ambient temperature 
Ta outside the cell, the electrolyte will exchange heat with the 
environment surrounding it through the walls of the container in 
which it is stored. The thermal power exchanged T Sext will be 
given by: 


T Sexi = HA(T, — T) (23) 


where H is the overall heat transfer coefficient and A is the con- 
tainer’s external surface area. If the thermal resistance due to 
film coefficients of the container’s external and internal walls is 
neglected, H will stand for the heat conduction transfer: 


H=- (24) 
E 
where k is the thermal conductivity of the material and e the 
thickness of the container walls. 
Fig. 9 shows the complete block diagram of the battery taking 
account of the above effects. 


4. Simulation results 


The model developed has been used to represent several situa- 
tions in order to observe their behavior, particularly with regard 
to thermal effects, as well as, for contrasting the results with 
some of those found in literature. Annex I contains the values 
taken into account in the simulation as well as the references 
from which they were taken. 

Firstly, the charge process of a six element cell coupled in 
series (corresponding to a standard 12 V. battery) was simulated 
for two different charge voltages. 

Figs. 10-12 show the results of charging a battery at a volt- 
age of 14 V for 2500 s, with an initial electrolyte temperature of 
298.15 K. The process was applied to an isothermal model and 


=n 


Leh, 

$) 

: 2 eee |. T 
P](9+12)+(20)f ] 

___Si" “Cathode main reaction $ 

———— a (+) 


i E(0) a 
E (5)+(19) PIi 
A E(i) 
T 
k a 


| _|(9+12)+(20)] | 


Fig. 9. Complete block diagram of the battery. 
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Fig. 10. Temperatures for a 14 V charge. Thick line, adiabatic thermal model. 
Thin line, non-adiabatic thermal model. 
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Fig. 11. State of Charge for 14 V charge. Thick solid line, isothermal model. 
Thin solid line, non-adiabatic thermal model. Dashed line, adiabatic thermal 
model. 
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Fig. 12. Charge currents for 14 V charge. Thick solid line, isothermal model. 
Thin solid line, thermal model. Dashed line, adiabatic thermal model. 
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Fig. 13. Temperatures for a 12.2 V charge. Thick solid line, adiabatic thermal 
model. Thin solid line, non-adiabatic thermal model. 


two thermal models. In a thermal model, which we have called 
the adiabatic model, the battery does not interchange thermal 
energy with the outside. In the other thermal model, the battery 
interchanges thermal energy with the outside and we have called 
it the non-adiabatic model. It can be seen that there are consid- 
erable differences depending on the model chosen. In the case 
where thermal exchange with the outside was taken into con- 
sideration, the electrolyte temperature underwent an increase in 
temperature of 11 K, and 21 K in the other case. The state of 
charge (SOC) values vary between 94 and 84% at the end of the 
charge process for both extreme cases. The charge current value 
at the end of the process is 55 A for the isothermal model, 35.8% 
greater than for the adiabatic model and 15.1% greater than for 
the non-adiabatic model. 

Likewise, Figs. 13—15 show the results obtained during charg- 
ing at a voltage of 12.2 V. As can be seen, in this case the charge 
current is much lower, and as a result the differences among the 
three models are therefore negligible. The difference found in 
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0 2000 2500 


500 1000 1500 
Time (sec.) 
Fig. 14. State of charge for a 12.2 V charge. Thick solid line, isothermal model. 


Thin solid line, non-adiabatic thermal model. Dashed line, adiabatic thermal 
model. 
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Fig. 15. Charge currents for a 12.2 V charge. Thick solid line, isothermal model. 
Thin solid line, non-adiabatic thermal model. Dashed line, adiabatic thermal 
model. 


the charge current values at the end of the process was no greater 
than 3% in the worst case. 


5. Model validation 


A model of these characteristics is highly dependent on a large 
number of parameters, not only of an electrochemical nature, 
but also of a constructive nature (initial electrolyte volume and 
density, size and internal surface area of the electrode, thermal 
characteristics of the container material, its surface area in con- 
tact with the air, its thickness, . . .). It is, therefore, complicated 
to validate the model from either experimental or theoretical 
results found in literature, as the cell or model used as a refer- 
ence would need to be perfectly defined. It is almost impossible 
to do this by taking published data. A correct verification of the 
model would, therefore, require a series of laboratory tests to be 
carried out. This will be the aim of a future work. 

Nevertheless, in order to confirm the good behavior of the 
model it is necessary that some comparison should be brought to 
bear with either experimental results or those provided by other 
sufficiently contrasted models. Ceraolo [7] proposes a model for 
lead-acid batteries, developed by constructing equivalent elec- 
tric circuits. This model has been validated with laboratory tests 
provided by the author himself [8]. This model and the data 
used for its validation have been used to transfer them to the 
model developed in this paper and be able to proceed to an 
initial validation. 

The battery used in [7] has a Cio capacity of 540 Ah. Once 
completely charged and stabilized, it was discharged by a 63 A 
current for a period of 7.197 h, at the end of which time the cur- 
rent remained zero. In order to validate the model developed, 
the same process was reproduced using the non-adiabatic ther- 
mal model presented in this paper. The results relative to the 
voltage per element supplied by the battery during the discharge 
process under the conditions described are shown in Fig. 16. As 
can be seen from Table 1, they coincide fairly closely with those 
supplied by the experimental results obtained in Ref. [8]. 


2.15 
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Fig. 16. Discharge voltage for a 63 A current, followed by a period of zero 
current. 


Table 1 
Voltage values during the discharge process (see Fig. 10) 

Vo (V) Vı (V) V2 (V) V3 (V) 
Ref. [8] 2.115 1.995 2.01 1.788 
This paper 2.119 1.962 2.02 1.786 


A discharge at constant current was simulated up to a 1.8 V 
cut-off voltage per element, in periods of 1, 3 and 10h (Fig. 17), 
in order to obtain the values for the capacities for C1, C3 and 
C10 and to be able to compare them with those provided by the 
Ref. [8]. As can be seen in Table 2 and in Fig. 18, the values 
coincide fairly well. 
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Fig. 17. Discharge voltage at constant current for J4, /3 and Jo. 


Table 2 
Data relative to capacity for different discharge rates 

Cı (Ah) C3 (Ah) Cio (Ah) 
Ref. [8] 227.5 357 540 
This paper 233 361.5 457 
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Fig. 18. Capacity vs. discharge time. 
6. Conclusions 


This work has presented a model of an electrochemical cell 
including complex behavior, in order to show the different 
behavior presented in the charging—discharging process of an 
electrochemical cell. 

The charging process for a standard 12 V battery for two dif- 
ferent voltage charges was simulated using both the isothermal 
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model and the thermal model. With the latter, it can be seen that 
when the charge currents are high, the thermal effects are con- 
siderable, and at the end of the charge process the electrolyte 
reaches temperatures considerably higher than at the beginning. 
This increase in electrolyte temperature has repercussions on 
the charge current value, and as a result, on the state of charge. 
Important mistakes can be made if an isothermal model is used. 
Also, important differences can be seen according to whether the 
cell’s heat exchange with the environment is taken into account 
or not. It is, therefore, important to make an efficient represen- 
tation of the heat transmission process between the electrolyte 
and the environment surrounding the cell. 

For those processes where the charge currents are moder- 
ate, the thermal effects are not significant, and the differences 
between both models negligible. This means that in many sit- 
uations a much simpler isothermal model can be used without 
committing any significant errors. 
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Parameter Value Ref. 
Initial specific surface area for Pb electrode amax,Pb (cm7!) 23,000 [20] 
Initial specific surface area for PbO? electrode amax,PbO, (cm7!) 230,000 [20] 
Initial exchange current density io,max;Pb Eq. (2) (A cm~?) 4.96 x 106 [20] 
Initial exchange current density imax,Pbo, Eq. (3) (A cm~?) 3.19 x 1077 [20] 
Initial exchange current density imax,H, Eq. (14) (A/cm) 6.607 x 107!4 [24] 
Initial exchange current density imax,o, Eq. (17) (A cm~?) 1.3 x 10714 [16] 
Initial cell temperature T (K) 298.15 — 
Ambient temperature T, (K) 298.15 — 
Standard reduction potential Ep, Eq. (2) (V) —0.356 [25] 
Standard reduction potential Ep, Eq. (3) (V) 1.685 [25] 
Standard reduction potential En, Eq. (14) (V) 0.000 [25] 
Standard reduction potential Eò, Eq. (17) (V) 1.229 [25] 
Charge transfer coefficient app Eq. (2) 0.775 [20] 
Charge transfer coefficient «ppo, Eq. (3) 0.575 [20] 
Charge transfer coefficient ay, Eq. (14) 0.42 [24] 
Charge transfer coefficient «o, Eq. (17) 0.657 [16] 
Exponent for area correction’, Pb electrode ôpp 1.5 [5] 
Exponent for area correction*, PbO? electrode dppo, 1.5 [5] 
Exponent for exchange current density correction? Epp Eq. (2) 0.0 [17] 
Exponent for exchange current density correction? &ppo, Eq. (3) 0.3 [17] 
Exponent for exchange current density correction? ém Eq. (14) 1.0 [17] 
Exponent for exchange current density correction? o, Eq. (17) 0.0 HA 
Pb electrode size™® Vpp (cm?) 49 - 
Pb electrode size®-f VPbo, (cm3) 151.9 — 
PbO; electrode size™® Vy, (cm?) 49 — 
PbO, electrode size”! Vo, (cm?) 151.9 — 
Mean ionic activity coefficient y+ Fitted from ref. data [22] 
Water activity XH,0YH0 Fitted from ref. data [22] 
Specific heat capacity of electrolyte cy (J Kg7! °C7!) Fitted from ref. data [22] 
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Parameter Value Ref. 
Temperature coefficient? dE},/0T Eq. (2) (from $° tables) (V K7!) —1.154 x 1073 [25] 
Temperature coefficient? 3Eppo, /3T Eq. (3) (from S° tables) (V K!) 3.27 x 1074 [25] 
Temperature coefficient? 3E, /aT Eq. (17) (from $° tables) (V K7!) —8.46 x 1074 [25] 
Thermal conductivity of container k (Wm-!°c"!) 0.2 — 
Thickness of container walls ¢ (m) 0.005 = 
External surface area of container? A (m?) 0.26 — 
External surface area of container? A (m?) 1.03 = 
Electrolyte volume/element® (1) 1.0 - 
Electrolyte volume/element! (1) 5.25 = 
Maximum electrolyte density®f (Kg17!) 1.28 = 
Maximum electrolyte concentration [SO4H? |max (mol Kg"!) 6.08 — 


* The specific surface area for each electrode changes with SOC so that: aj = dmax,;(1 — Q/Qmax)ôj for charge, and aj = dmax, ;(Q/Qmax)ôj for discharge; 


j= Pb,PbO2. 


b The exchange current density changes for each reaction with the concentration of acid [SO4H2], so that: io, ; = io,max, ;([SO4H2]/ [SO4H2 Imax) j= 


Pb, PbO2, Ho, Oo. 


€ Current intensity / for any reaction is related to the corresponding current density i, the specific surface area, and the size of the electrode V by means of: J; = ija; Vj; 


J=Pb, PbOz, H2, O2. 


d To calculate Srey using expression (19), taking as reference value that corresponding to Eq. (14). 


e Values corresponding to the results in Figs. 10-15. 
f Values corresponding to the results in Figs. 16-18 and Tables 1 and 2. 
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